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Abstract: Adiabatic quantum computation for performing quantum computations such as 
Shor's algorithm is protected against thermal errors by an energy gap of size 0(l/n), where 
n is the length of the computation to be performed. 



Adiabatic quantum computing is a novel form of quantum information processing that 
allows one to find solutions to NP-hard problems with at least a square root speed-up [1]. 
In some cases, adiabatic quantum computing may afford an exponential speed-up over 
classical computation. It is known that adiabatic quantum computing is no stronger than 
conventional quantum computing, since a quantum computer can be used to simulate an 
adiabatic quantum computer. Aharonov et al. showed that adiabatic quantum computing 
is no weaker than conventional quantum computation [2]. This paper presents novel models 
for adiabatic quantum computation and shows that adiabatic quantum computation is 
intrinsically protected against thermal noise from the environment. Indeed, thermal noise 
can actually be used to 'help' an adiabatic quantum computation along. 

A simple way to do adiabatic versions of 'conventional' quantum computing is to use 
the Feynman pointer consisting of a line of qubits [3]. The Feynman Hamiltonian is 
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H = - U e ® \£ + 1)(£\ + U\ ® \£)(£ + 1|, (1) 
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where Ut is the unitary operator for the £'th gate and \£) is a state of the pointer where 
the tth qubit is 1 and all the remaining qubits are 0. Clearly, H is local and each of its 
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terms acts on four qubits at once for two-qubit gates. If we consider the pointer to be a 
'unary' variable, then the each of the terms of H acts on two qubits and the unary pointer 
variable. 

Assume that the computation has been set up so that all qubits start in the state 
0. The computation takes place and the answer is placed in an answer register. Now a 
long set of steps, say n/2 takes place in which nothing happens. Then the computation is 
undone, returning the system to its initial state at the n — l'th step. The computational 
circuit then wraps around and begins again. The eigenstates of H then have the form 

n-l 

\b, k) = (1/y/n) e t27rhi/n U e ...U \b)® \£), (2) 

with eigenvalue -2 cos 27r/c/n, since H\b,k) = -( e -« 27rA: /"- _|_ e l2 ' Kk / n )\b, k). These can be 
thought of as momentum eigenstates for the propagation of the pointer qubit down the 
chain. The |0, k) momentum eigenstates have the nice feature that if you measure the 
answer register, the probability of obtaining the answer is 1/2. 

Feynman used this Hamiltonian to set up a traveling state of the pointer (a Gaussian 
superposition of the momentum eigenstates) so that the computation could take place 
sequentially. This propagating state is mathematically equivalent to a coherent quantum 
walk down the chain [4-5]. Landauer pointed out that for any realistic implementation of 
such a system, imperfections in the Hamiltonian would result in Anderson localization and 
the computation would get 'stuck' and fail to propagate to the end [6]. Localization is a 
significant problem for quantum walks in general and for the Feynman quantum computing 
model in particular. 

One can also use Feynman's Hamiltonian to implement the computation adiabatically. 
In this will be seen, localization does not affect the computation. Let 

if = ^(l-|0) J (0|)®|£ = 0)(£ = 0| (3) 

3 

be the Hamiltonian one of whose ground states is the state with all the computer qubits 
equal to zero and the pointer at the zero spot of the line. Hq is degenerate: states 
where the pointer is at other places also have zero energy, regardless of the values of the 
computational qubits. First, prepare the system in the state |00. . .0)\£ = 0) by applying 
the Hamiltonian r]H + H 1: where Hi = —\i = 0)(£ = 0| and 1] is a positive term that sets 
the overall energy scale. 
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Now gradually turn on an H term while turning off the Hi term: the total Hamiltonian 
is 7]H + (1 — X)H 1 + XH. As A is turned on, the system goes to its new ground state 
\b = 0, k = 0). It can be verified numerically and analytically [7] that the minimum energy 
gap in this system occurs at A = 1: consequently, the minimum gap goes as 1/n 2 . In fact, 
the energy gap due to the interaction between the Hi and the H terms is just the energy 
gap of the simpler system consisting just of the chain qubits on their own, confined to the 
subspace in which exactly one qubit is 1: that is, it is the energy gap of a qubit chain 
with Hamiltonian (1-X)H 1 + XH > , where H' = -^2 e \e+l)(£\ + \f)(e+l\. This gap goes 
as 1/n 2 . Accordingly, the amount of time required to perform the adiabatic passage is 
polynomial in n. 

When the adiabatic passage is complete, the energy gap of the H term on its own 
goes as 1/n 2 from the cosine dependence of the eigenvalues of H: it is also just the 
energy gap of the simplified system in the previous paragraph. This implies that the 
adiabatic passage can accurately be performed in a time polynomial in n. Measuring the 
answer register now gives the answer to the computation with probability 1/2. This is 
an alternative (and considerably simpler) demonstration to that of [2] that 'conventional' 
quantum computation can be performed efficiently in an adiabatic setting. 

An interesting feature of this procedure is that the adiabatic passage can be faulty 
and still work just fine: all of the energy eigenstates in the \b = 0, k) sector give the correct 
answer to the computation with probability 1/2, for any k. The real issue is making sure 
we do not transition to the \b ^ 0,k) sector. But the Hamiltonians Hi and H do not 
couple to this sector: so in fact, we can perform the passage non-adiabatically and still 
get the answer to the computation. For example, if we turn off the Hi Hamiltonian very 
rapidly and turn on the H Hamiltonian at the same time, the system is now in an equal 
superposition of all the \b = 0, k) eigenstates. If we wait for a time oc n 2 (corresponding 
to the inverse of the minimum separation between the eigenvalues of H) , then the state of 
the system will be spread throughout the \b = 0, k) sector, and we can read out the answer 
with probability 1/2. This method effectively relies on dispersion of the wavepacket to find 
the answer. 

Since coherence of the pointer doesn't matter, we can also apply a Hamiltonian to 
the pointer that tilts the energy landscape so that higher pointer values have lower energy. 
(E.g., -^pointer = — ^£ £E\£) Starting the pointer off in the initial state above and 
letting it interact with a thermal environment will obtain the answer to the computation 
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in time of 0(n). Similarly, in the absence of an environment, starting the pointer off in a 
wavepacket with zero momentum at time and letting it accelerate will get the answer to 
the computation in even shorter time. 

Clearly, this method is quite a robust way of performing quantum computation. Let 
us look more closely at the sources of this robustness. If rj is big, then there is a separation 
of energy of 0(r]/n) between the \b = 0, k) sector — states which give the correct answer 
to the computation — and the \b ^ 0, k) sector — states which give the incorrect answer 
to the computation. This is because {b ^ 0, k\i]Ho\b ^ 0, k) = rj/n. This energy gap goes 
down only linearly in the length of the computation and can be made much larger than 
the gap between the ground and first excited state by increasing r\ » 1. 

This second energy gap is very useful: it means that thermal excitations with an 
energy below the gap will not interfere with obtaining the proper answer. That is, this 
method is intrinsically error-tolerant in the face of thermal noise. Indeed, it is this 0(i]/n) 
gap that determines how rapidly the computation can take place rather than the 0(l/n 2 ) 
gap between the ground and excited states. 

Of course, the actual errors in a system that realizes the above scheme are likely to arise 
from variability in the applied Hamiltonians. The energy gap arguments for robustness 
only apply to the translational dynamics of the system (this is what makes the analysis of 
the system tractable in the first place) . That is, errors that affect each Ug on its own are 
not protected against: but these are the errors that cause the computation to come out 
wrong. Of course, one can always program the circuit to perform 'conventional' robust 
quantum computation to protect against such errors. One must be careful, however, that 
errors that entangle the pointer with the logical qubits do not contaminate other qubits: 
conventional robust quantum computation protocols will have to be modified to address 
this issue. Farhi et al. have recently exhibited error correcting codes for 'conventional' 
adiabatic quantum computation [1] that can protect against such computational errors [8] . 

The use of error correcting codes to correct the variation in the Ug may well be overkill. 
In any system manufactured to implement adiabatic quantum computing, these errors in 
the Ug are essentially deterministic: the Ug could in principle be measured and their varia- 
tion from their nominal values compensated for by tuning the local Hamiltonians. Because 
it involves no added redundancy, such an approach is potentially more efficient than the 
use of quantum error correcting codes. Exactly how to detect and correct variations in the 
Ug will depend on the techniques (e.g., quantum dots or superconducting systems) used to 
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construct adiabatic quantum circuits. 

It is also interesting to note that performing quantum computation adiabatically is 
intrinsically more energy efficient than performing a sequence of quantum logic gates via 
the application of a series of external pulses. The external pulses must be accurately 
clocked and shaped, which requires large amounts of energy. In the schemes investigated 
here, the internal dynamics of the computer insure that quantum logic operations are 
performed in the proper order, so no clocking or external pulses need be applied. The 
adiabatic technique also avoids the Anderson localization problem raised by Landauer. 

The above construction requires an external pointer and four qubit interactions. One 
can also set up a pointerless model that requires only pairwise interactions between spin 1/2 
particles (compare the following method with the method proposed in reference [9]). Let 
each qubit in the computational circuit correspond to a particle with two internal states. 
Let each wire in the circuit correspond to a mode that can be occupied by a particle. The 
€th quantum logic gate then corresponds to an operator Hg = Ag + A\, where Ag is an 
operator that takes two particles from the two input modes and moves them to the output 
modes while performing a quantum logic operation on their two qubits. That is, 

A i = al utl a inl al ut2 a in2 <E> Ug. (5) 

Note that Ag acts only when both input modes are occupied by a qubit-carrying particle. 
If we use the Hamiltonian H = J2e Hg in place of H in the construction above, the 
ground state of this Hamiltonian is a superposition of states in which the computation 
is at various stages of completion. Just as above, measurement on the ground state will 
reveal the answer to the computation with probability 1/2. 

Note that even though the Hamiltonian in equation (5) involves a product of opera- 
tors on four degrees of freedom (the internal degrees of freedom of the particles together 
with their positions), it is nonetheless a physically reasonable local Hamiltonian involv- 
ing pairwise interactions between spin-1/2 particles. To simulate its operation using an 
array of qubits as in [2] would require four qubit interactions, as in the pointer model dis- 
cussed above. This point is raised here because of the emphasis in the quantum computing 
literature on reducing computation to pairwise interactions between qubits. Pairwise in- 
teractions between particles or fields - i.e., the sort of interactions found in nature - may 
correspond to interactions between more than two qubits. 

Without further restrictions on the form of the quantum logic circuit, evaluating the 
energy gap in this particle model is difficult, even for the final Hamiltonian 1 — H. But we 
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can always set up the computational circuit in a way that allows the adiabatic passage to 
be mapped directly on the Feynman pointer model above. The method is straightforward: 
order the quantum logic gates as above. Now insert additional quantum logic gates between 
each consecutive pair of gates in the original circuit. The additional gate inserted between 
the tth and £+Vth quantum logic gates couples one output qubit of the £'th quantum logic 
gate with one input qubit of the t + l'th gate, and performs a trivial operation U = 1 on 
the internal qubits of these gates. The purpose of these gates is ensure that the quantum 
logic operations are performed in the proper sequence. Effectively, one of the qubits from 
the tth gate must 'tag' one of the qubits from the I + l'th gate before the I + l'th gate 
can be implemented. Accordingly, we call this trick, a 'tag-team' quantum circuit. 

Tag-team quantum circuits are unitarily equivalent to the Feynman pointer model 
with an extra, identity gate inserted between each of the original quantum logic gates. 
Accordingly, the spectral gap for tag-team quantum circuits goes as 1/n 2 and the quantum 
computation can be performed in time 0(poly(n)). Just as for the pointer version of 
adiabatic quantum computing, the important spectral gap for tag-team adiabatic quantum 
computation is not the minimum gap, but rather the gap of size 0(rj/n) between the 
ground-state manifold of 'correct' states and the next higher manifold of 'incorrect' states. 
Once again, the existence of this gap is a powerful protection against errors in adiabatic 
quantum computation. 

The methods described above represent an alternative derivation of the fact that 
adiabatic quantum information processing can efficiently perform conventional quantum 
computation. The relative simplicity of the derivation from the original Feymnan Hamil- 
tonian [3] allows an analysis of the robustness of the scheme against thermal excitations. 
Adiabatic implementations of quantum computation are robust against thermal noise at 
temperatures below the appropriate energy gap. The appropriate energy gap is not the 
minimum gap, which scales as n 2 , but the gap between the lowest sector of eigenstates, 
which give the correct answer, and the next sector. This gap scales as rj/n, where n is an 
energy parameter that is within the control of the experimentalist. 
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